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ABSTRACT 

I present three methods to determine the distance to the Galactic centre Rq, the solar 
azimuthal velocity in the Galactic rest frame V 9} q and hence the local circular speed 
V c at Rq. These simple, model-independent strategies reduce the set of assumptions 
to near axisymmetry of the disc and arc designed for kinematically hot stars, which 
are less affected by spiral arms and other effects. 

The first two methods use the position-dependent rotational streaming in the 
heliocentric radial velocities (U). The resulting rotation estimate 9 from U velocities 
does not depend on V^q. The first approach compares this with rotation from the 
galactic azimuthal velocities to constrain V^© at an assumed Rq. Both V 3i q and Rq 
can be determined using the proper motion of Sgr A* as a second constraint. The 
second strategy makes use of 9 being roughly proportional to Rq. Therefore a wrong 
Rq can be detected by an unphysical trend of V 9: q with the intrinsic rotation of 
different populations. From these two strategies I estimate Rq = (8.27 ± 0.29) kpc 
and V 9t Q = (250 ± 9)kms _1 for a stellar sample from SEGUE, or respectively V c — 
(238 ± 9)kms~ 1 . The result is consistent with the third estimator, where I use the 
angle of the mean motion of stars, which should follow the geometry of the Galactic 
disc. This method also gives the Solar radial motion with high accuracy. 

The rotation effect on U velocities must not be neglected when measuring the Solar 
radial velocity Uq. ft biases Uq in any extended sample that is lop-sided in position 
angle a by of order lOkrns -1 . Combining different methods I find Uq ~ I4kms _1 , 
moderately higher than previous results from the Geneva-Copenhagen Survey. 
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1 INTRODUCTION 

Among the central questions in Galactic structure and pa- 
rameters are the Solar motion, Solar Galactocentric radius 
Ro and the local circular velocity V c of our Galaxy. Galac- 
tic rotation curves are found to be gen erally quite flat over 
a va st range of Galactocentric radii |Krumm fc Salpeterl 
1 19791). As common for Galaxies with exponential discs 



1 Freemadfl970h there is some evidence for a radial trend of 
the Galacti c circular velocity near th e Sun, but it is very 
mode rate (|Feast fc Whitelockl Il997l ; iMcMillan fc Binnevl 
2010), so that the circular velocity at solar Galactocentric 
radius Ro characterises well the entire potential. 

Initially local kinematic data from stellar samples were 
the main source to extract Galactic parameters including 
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V c , e.g. by the use of the Oort constants (|Oort| [l927'). The 
kinematic heat of stellar populations requires large sample 
sizes, so that studies determining the Local Standard of Rest 
(LSR) are still primarily based on stars, but some classi- 
cal strategies like th e position deter mination of the Galactic 
centre pioneered by IShaplevI (|l918h reached their limits by 
the constraints from geometric parallaxes, by the number of 
available luminous standard candles, and by the magnitude 
requirements of stellar spectroscopy. 

Apart from some more recent att empts to use lumi- 
nous stars (e.g. iBurton fc Ba nia 1974), most current ev- 
idence on Ro and V c de rives from mode l ling streams in 
the G alactic halo (see e.g. Ilbata el al.ir200ll ;fMaicwsk i~et al.l 
2006) and from radio observations of the HI terminal ve- 
locity (see lMcMillanll201ll h the Galactic centre, m olecular 
clouds and MASERs (e.g. iReid fc Brunthaleill2004h . While 
the first branch relies strongly on assumptions such as dis- 
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tance scale and shape of the Gal actic potential (cf. the dis- 
cussion in iMaiewski et alj|2006l ). there are further compli- 
cations hte_^i«_^fhrreof_tidal streams to delineate stellar 
orbits (|Evre fc Binnevll2009h ■ 

Radio observations l|Reid fc Brunth aler 2004) have ac- 
curately determined the proper motion of the radio source 
Sgr A*, which is identified with the central black hole of 
the Milky Way (for a discu ssion of possible uncertainties see 
also iBroderick et alj|20ilh . It tightly constrains the ratio 
of the solar speed in the azimuthal direction V 9 ,0 to Ro, 
but further information is required to obtain both quanti- 
ties, letting aside the need for an independent measurement. 
Recently parallaxes to objects in the central Galactic re- 
gions have become available (see the discussion in lReid et al.l 
I2009T ) . and values for the Galactic circular spee d have been 
deriv ed from HI motions and from MASERs (Rvgl et al. 

Despite the decent errors in the determined kinemat- 
ics of MASERs, the small sample sizes impose considerable 
systematic uncertainties: they are not on pure circular orbits 
and more importantly they are intimately connected to the 
intense star formation in spiral arms, where the kinematic 
distortions are largest. As the youngest strategy, studies of 
orbits in the Galactic centre have gained high precision, but 
weakly constrai n Rq due to a str ong degeneracy with the 
black hole mass (|Ghez et alj|2009h . Hence the values for Rq 
and V g ,o remain under debate. 

An independent determination of Galactic parameters 
is facilit ated by the new larg e spectroscopic surveys like 
RAV E (|Steinmetz et all 120051 ) and SEGUE (|Yannv et al.l 
2009). I will show that already now the stellar samples, 
which so far h ave been primarily used for the explorat ion of 
substructure (jBelokurov et al]|2007l ; Ifiahn et all 1201 lh give 
results for V c competitive with radio observations. 

Common ways to determine Galactic parameters from 
the motion of stars require a significant bundle of assump- 
tions and modelling efforts to evaluate the asymmetric drift 
in a subpopulation compared to the velocity dispersion and 
other measurements and assumptions, e.g. detailed angular 
momentum and energy distributions, the validity of the the- 
oretical approximations, etc. On the contrary a good mea- 
surement strategy for Galactic parameters should have the 
following properties: 

• Do not rely on dispersions and other quantities that 
require accurate knowledge of measurement errors. 

• Do not rely on kinematically cold objects in the disc 
plane, which are particularly prone to perturbations from 
Galactic structure. 

• Do not rely on any specific models with hidden assump- 
tions and parameters and require as few assumptions as pos- 
sible. 

In an attempt to approximate these conditions this note 
concentrates on the use of kinematically hot (thick disc, 
halo) populations. The use of mean motions and the as- 
sumption of approximate axisymmetry in our Galaxy will 
be sufficient to fix the Galactocentric radius of the Sun as 
well as the Solar motion and Galactic circular speed. 

In Section[2]I lay out the general method before describ- 
ing the sample selection and treatment in Section[3] The lat- 
ter comprises a discussion of proper motion systematics, a 
discussion of distances in subsection [3]3] and a discussion of 
line-of-sight velocities and the vertical motion of the Sun. In 
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Figure 1. Definition of kinematic and geometric quantities. On 
the connection line between Sun and Galactic centre (abbrevi- 
ated GC) the heliocentric (U, V, W) and galactocentric velocities 
(Ug, V g , Wg), i.e. velocities in the galactocentric cylindrical coor- 
dinate system are identical, while in general they differ. This gives 
rise to systematic streaming motion in the heliocentric frame. 



Section [4] I lay out the radial velocity based rotation mea- 
surement using SEGUE, discuss the radial velocity of the 
Sun in Subsection 14.11 and shortly discuss possible trends. 
From 14.31 on I present three methods to to extract Galac- 
tic parameters from stellar samples, the first two methods 
relying on the radial velocity based rotation measurement 
and the third on the mean direction of stellar motions. In 
Sectional summarise the results. 



2 GENERAL IDEA 

Before I start the dirty work of data analysis it seems appro- 
priate to concisely lay out the general definitions and ideas. 



2.1 Definitions 

This work is based on comparisons between velocities in the 
heliocentric frame and galactocentric velocities. In the helio- 
centric frame I define the velocities U, V, W in a Cartesian 
coordinate system as the components pointing at the Sun 
towards the Galactic centre, in the direction of rotation and 
vertically towards the Galactic North pole. The galactocen- 
tric velocities are defined in cylindrical coordinates around 
the Galactic centre with U g ,V g , W g pointing again towards 
the Galactic central axis, in the direction of rotation and ver- 
tically out of the plane, as illustrated in Fig. [l] Uq, V®, Wq 
are the three components of the Solar motion relative to the 
Local Standard of Rest (i.e. the circular orbit at the local 
galactocentric radius R p), for which I assume the values from 
ISchonrich etafl i|2010l ) if not stated otherwise. The total az- 
imuthal velocity of the Sun in the Galactic frame is written 
Vg,© = Vq + Vc where V c denotes the circular velocity in the 
disc plane of the Milky Way at Rq. The mean rotation speed 
of a subsample of stars will be named 8 and by the asymmet- 
ric drift is generally smaller than V c - For convenience local 
Cartesian coordinates are defined left-handed x,y,z point- 
ing outwards, in the direction of rotation and upwards with 
the Sun at the origin. 
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2.2 An absolute measure of rotation 

Commonly studies on the rotation of Galactic components 
suffer from the uncertainty in the solar azimuthal velocity 
that directly translates into a systematic uncertainty of stel- 
lar azimu thal velocities. T h ere is a way out: As already dis- 
cussed in ISchonrich et alJ (|201ll ) and evident from Fig. [T] 
Galactic rotation leaves its imprint in both the heliocentric 
U and V velocities. The direction of the rotational compo- 
nent of motion V g turns with the angle a between the lines 
Galactic centre - Sun and Galactic centre - star to align 
partly with the radial component U of the local Cartesian 
frame. From this systematic streaming motion in heliocen- 
tric U velocities we can directly infer the Galactic rotation of 
a stellar sample. Accounting for rotation the mean motions 
of stars in the heliocentric frame are: 



Sun' :Sun 



(U) = -U Q +6sma 

(^} = -V s , Q +0(l-cosa) 

(W) = -W e , 



(1) 



with the rotation speed of the population 9 and the Galactic 
angle a. (Uq, V g ,Q, Wq) are the velocity components of the 
Solar motion in the Galactic rest frame radially towards the 
Galactic centre, azimuthally in the direction of disc rotation 
and vertically out of the plane. Observing the change of he- 
liocentric U velocity in a sufficiently extended sample hence 
estimates the rotation of a component 8, once we know the 
angle a. 

Similarly the observed heliocentric azimuthal velocity 
V shrinks towards larger |sinaj. Prima facie it might be 
tempting to use this term. Yet, all it expresses is a slowing 
of heliocentric azimuthal velocities. In general 8 is not a real 
constant, but a function of altitude z, galactocentric radius 
R, metallicities, etc. In particular the kinematically hotter 
disc populations at higher altitudes above the p lane have 
a larger asymmetric drift (see e.g. iBinnevi [2010h . or vice 
versa slower rotation, seconded by an increase in the share 
of the nearly non-rotating halo. Pointing out from the Sun's 
position in the mid-plane, | cosa correlates via the typical 
distance range with altitude z. In contrast to U, where this 
bias is of second order, it is of first order in V and prohibits 
any naive estimate of 8 from azimuthal velocities. A similar 
uncertainty derives from the possibility for variations of 8 in 
R even at the same altitude. 

The radial velocities are better behaved. Sample 
changes along the baseline sin a are not only a second order 
effect, but are generally more robust: the greatest concern is 
a bias in the derived Uq, which can be prevented by fixing 
its value by previous knowledge or by simply having a suf- 
ficient sky coverage that allows to balance both sides of the 
Galactic centre and includes sufficient numbers of stars at 
low | sin«|. If the integrity of Uq is ensured, sample changes 
along the baseline just result in different statistical weights 
of the involved populations. 

The value 8 derived from heliocentric radial motions 
is by construction independent from the assumed Solar az- 
imuthal velocity and (in an extended sample) is only weakly 
impacted by (if held fixed) a mis-estimate of Uq. This con- 
trasts to direct measurements of V g that rely fully on the 
total solar azimuthal velocity V g ,Q - a short-come now turn- 
ing into a virtue: V g ,® and Ro can be varied to force an 
agreement between the rotation estimates from 8 and V g . 




Figure 2. Changing the Galactocentric radius of the Sun Ro. 
The Galactic longitude of the star does not change so that the 
angle a grows with shrinking Rq . This reduces the rotation speed 
estimated from heliocentric radial velocities proportionally. 



As the position-dependent statistical weights in 8 are known, 
one can estimate a weighted average of V g exactly on these 
weights. The sole assumption of approximate axisymmetry 
of the Galactic disc (but not necessarily in the sample se- 
lection) is then sufficient to compare 8 and V g and measure 
Galactic parameters without the least concern about sam- 
ple selection, velocity trends, etc. The comparison gives a 
relation between V 9: q and Tio, which can be combined with 
independent data like the proper motion of Sgr A* to pro- 
vide a superior estimate of Galactic parameters. However, 
Ro alone can be measured with this simple method alone, as 
6 is roughly proportional to Ro (see Section [4.31 or Fig. [2}, 
so that any error in Tio shears the estimate of V g ,Q for fast 
rotating populations against their slower counterparts. 



3 SAMPLE SELECTION AND DISTANCES 
3.1 Data selection 

Any reader not interested into the details of sample selec- 
tion, distance analysis and the proper motions may skip this 
section and turn directly towards Section [4. II 

All data used in this study were taken from the sev- 
enth and eighth da t a release (hereafte r DR7 and DR8 
lAbazaian et al.ll2009l ; lAihara et al"1l201ll) of th e Sloan Dig- 
ital Sky Survey fSDSS lEisenstein et all l201lh . The stellar 
spectra stem from the Sloan Extens ion for Galactic U nder- 
standing and Exploration (SEGUE. jYannv et al.ll2009l '). An 
evaluation of the perfo rmance of their p arameter pipeline 
(SSPP) can be found in lLee et al] (|2008al lbh. This work has 
no interest in specific metallicity distributions, but in hav- 
ing kinematically unbiased samples that include as many 
stars as possible. Hence the sample drawn from DR8 con- 
sists of a raw dataset of 224 019 stars from all SEGUE tar- 
get categories and the photometric and reddening standard 
stars that do not include any proper motion cuts and have 
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clean photometry. In particular I use the categories of F 
turnoff/sub-dwarf, Low Metallicity, F/G, G dwarf, K dwarf, 
M sub-dwarf stars from SEGUE1, as well as MS turnoff, 
Low Metallicity and the reddening and photometric stan- 
dard stars from all samples. The categories were selected by 
the target bitmasks of the database. Some of these categories 
overlap with kinematically biased selections and the survey 
descriptions are not clear about a possibility of direct (cross- 
selection of stars) or indirect (non-selection of stars that are 
targeted in a biased category) bias0 I made sure that the 
presented results are unchanged when selecting stars by the 
unique target flags in the SSPP. Also checks on leaving out 
e.g. K dwarfs that might overlap with the biased K giants 
did not provide evidence for a notable bias by the SEGUE 
target selection. The older DR7 forms a subset of 182627 
stars. 

To ensure sufficient quality of the involved proper mo- 
tions, I select only stars with a match in the proper mo- 
tion identifications (match = 1), the presence of at least 
5 detections (nFit > 5) and a good position determina- 
tion with ora,cfd ec < 500 mas as a compromise between 
iMunn et alJ l|2004l ) and currently used selection schemes. I 
further require the formal errors on the proper motions to 
fulfil o" Mb ,(T ( j i < 4masyr -1 , removing the tail of uncertain 
proper motion determinations and require line-of-sight ve- 
locity errors below 20kms _1 . 

Sample homogeneity is not a concern in this study, but 
to avoid issues with reliability of the isochrones in colour 
regions with insufficient coverage to check them by statistics, 
I select < (g — i)o < 1, which tosses out a handful of stars 
especially on the red side. I tested that narrower colour cuts 
would have no significant effect on the presented results. 

For our purpose it is most important to have a good dis- 
crimination between dwarf stars and other categories and to 
dispose of as many evolved stars above the main sequence as 
possible. Hence only stars are used that have values for grav- 
ity and metallicity determined by t he pipeline. The s ample 
must be restricted to dwarf stars, as lSchonrich et all (|201lh 
showed that stars with intermediate gravities in the pipeline 
are an indeterminable blend of dwarf, subgiant and giant 
stars, which renders their distance scatter and residual bi- 
ases detrimen tal to any investigation of their kinematics. 

Following Schonric h et al.l (|201ll ) I adopt a sloping grav- 
ity cut generally tighter than the usually applied log(g) > 
4.0, demanding log(g) > 4.2 + 0.15[Fe/H] and setting a con- 
stant limit of log(g) > 3.9 for [Fe/H] < -2 and log(g) > 4.2 
for [Fe/H] > 0. It is important to point out that this se- 
lection could especially on the metal-poor end still be irre- 
sponsibly lenient for studies of detailed stellar kinematics, 
as it will comprise a serious contamination with misidenti- 
fied (sub-)giants and does not bias sufficiently against the 
brighter turn-off stars. Indeed the distance corrections de- 
scribed below react to a tighter gravity cut in the metal- 
poor regime by up to 10%. Relying only in mean velocities 
the preference in this work is leaned towards larger sample 
numbers and for assurance I checked that a tighter selection 
widens the error margins, but does not alter my results. 

Quite frequently studies on SEGUE use a cut in the 



signal to noise ratio S/N > 10 to ensure reasonable parame- 
ters. Since the determination of stellar parameters is difficult 
in those noisy spectra, the detection of a mild drift in the 
derived quantities towards higher S/N does not surprise, 
motivating a cut at S/N > 15 where results become more 
stable. 

De-reddening is done by subtraction of the reddening 
vectors given in the DR8 catalogue from the observed mag- 
nitudes. To limit uncertainties and avoid problems with com- 
plicated three dimensional behaviour of the reddening, I re- 
move all objects with an estimated g magnitude extinction 
A g > 0.75 mag. This strips all stars in the low latitude fields 
from the sample. I tested that the exclusion does not have a 
significant impact on my results apart from mildly increas- 
ing the derived formal errors, but maintain this cut for the 
sake of sample purity. I exclude all stars flagged for spectral 
anomalies apart from suspected carbon enhancement. The 
internal consistency of the distance estimator against closer 
distance bins starts eroding for s > 3.6 kpc. Hence I use this 
distance limit instead of the more common s < 4.0 kpc. 

3.2 Proper motions 

For a precise measurement of Galactic rotation I am re- 
liant on a decent control of systematic errors on proper 
motions. Unbiased noise is of no concern as proper mo- 
tion errors translate linearly into velocities (see equ. 20 in 
ISchonrich et alJl201ll ) and hence mean velocities (on which 
this paper relies) are only affected by the mean value of 
proper motion errors0 After a short description of proper 
motion systematics I will measure them in equatorial coor- 
dinates in Section 13.2.11 and develop terms that are more 
natural to a contamination problem and Galactic rotation 
in S ection 15321 

lAihara et al.l rt2011erratuml ') reported significant prob- 
lems with the DR8 astrometry implying also problems with 
proper motions. Fortunately the revised astrometry has be- 
come available during refereeing stage of this paper and I 
will exclusively use the new data. It should be noted that 
also DR7 had major problems, e.g. the entire rerun 648 was 
clearly erroneous/contaminated with large effects on kine- 
matics. Its 15% share of the total sample would bring down 
my later measurement of Galactic rotation by as much as 
14kms . Most problems with Sloan proper motions have 
been adressed in the revised astrometry, but the large im- 
pact calls for an investigation of residual biases. 

Two main sources of systematic terms can be identi- 
fied: astrometric "frame-dragging" and chromatic aberra- 
tions. Despite the efforts to clean up the Galaxy sample used 
for astrometry from residual contamination with Galactic 
stars, the sample may still be subject to some astrometric 
"frame-dragging" . Qualitatively those stars shift the frames 
in the direction of Galactic rotation or respectively their mo- 
tion relative to the Sun, underestimating the real motion in 
the sample. As an alternative explanation l|Bond et al.ll2009h 
traced the significant declination dependent net proper mo- 
tions on qua sars back to chromatic a berration by Earth's 
atmosphere l|Kaczmarczik et al.ll2009T ): the angle at which 



1 cf. the comments of Schlcsingcr ct al . (2011) at the analogous 
problem of the metallicity distribution. 



2 The knowledge of proper motion errors themselves plays only 
a subordinate role in the distance correction 
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fitting /j RA 


fitting fiDEC 


error 


/i,RA ' 10 3 


0.33 


0.18 


0.16 


/..DEC ■ 10 3 


-0.59 


7.88 


0.50 




0.115 


-0.205 


0.033 



Table 1. Parameters for systematic proper motion terms when 
fitting the linear equation 1 13. 2. Il l to detect systematic proper mo- 
tions on the quasar sample. Both declination (right column) and 
right ascension (left column) show significant systematics. Terms 
with > 2a significance are printed in bold letters. 



name 


DR8 


""DR8 


Di?8( 9 _ r ) >0 


fT DR8,( 9 -r)>0 


fi 


-0.0585 


0.0051 


-0.0615 


0.0055 


Cl 


-0.0005 


0.011 


0.001 


0.012 


ai 


0.00415 


0.00053 


0.00434 


0.00057 


fb 


-0.0589 


0.0054 


-0.0604 


0.0057 




-0.159 


0.013 


-0.159 


0.014 


"b 


0.00370 


0.00067 


0.00418 


0.00072 



Table 2. Fit parameters for systematic streaming in the quasar 
proper motions in DR8 with and without cut for (g — r) > as 
described in the text. 



the stellar light passes through our atmosphere strongly de- 
pends on the declination and this should give rise to colour- 
dependent offsets. Further, deformations in the telescope can 
introduce minor distortions of the focal plane. 



3.2.1 Measuring proper motion systematics in equatorial 
coordinates 

While the existence of these aberrati ons appears plausible, I 
tested proper motions on the cleaner [Schneider et alj (|201Ch 
quasar sample. I fitted a simple linear function in right as- 
cension and declination to the proper motions from the re- 
vised DR8 astronomy: 



RA , DEC 

— = Ji.RA • -j h /i.DEC • — ; h Ci 



(2) 



masyr _± deg deg 

where i — RA, DEC denotes the two possible proper mo- 
tion components in right ascension (RA) and declination 
(DEC), / are the coefficients, c are free fitting constants. 
Using the combined fit diminishes the impact from omitted 
variable bias that might arise from the particular sample 
geometry (the explaining variables are mildly correlated). 
The fit parameters are shown in Table [1] There is a clear 
and significant detection of systematic proper motions of 
order 0.2masyr -1 . I also checked the same fit on the older 
DR8 astrometry, where the derived terms were even slightly 
smaller. The large term on /dec, dec can be expected from 
chromatic aberration. However, the strong effect on the right 
ascension is more difficult to explain by atmospheric influ- 
ences, but may still come from specific distortions on the 
telescope or if the observations were made in a very par- 
ticular manner that correlates the airmass of a stellar ob- 
servation significantly to right ascension. The sample size of 
just about 31500 quasars limits the possibility to dissect into 
subgroups. However, cutting in observed colours there is no 
detectable difference between quasars at (g — r) < 0.1 and 
(g — r) > 0.3 implying that in the colour range interesting for 
the stellar sample the bias is quite consistent. Derived values 
go only adrift for {g — r) < 0, which bears no importance to 
the stellar sample. 

3.2.2 Proper motion systematics linked to rotation 

The strong deviations on the right ascension motivate tests 
for astrometric "frame" dragging by stellar contamination. 
In the following I develop a very crude analytic approxima- 
tion for a proper motion bias in the (I, fo)-plane as it could 
arise from Galactic streaming - note also that this toy is not 
used in any other way throughout the paper apart from the 
proper motion corrections: candidates for contamination of 
the galaxy sample should preferentially populate a certain 



magnitude range and also have preferred colours resembling 
the galaxies used. It hence seems appropriate to assume for 
simplification some dominant distance for stars contaminat- 
ing the sample, placed at s = 1 kpc. The outcomes will any- 
way not be critically affected by this distance (which can 
alter the shape of the systematic bias once the stars are suf- 
ficiently remote for the small angle approximation to break 
down). On this shell I assume the stars to rotate depen- 
dent on their altitude z around the Galactic centre with 
8(z) = (215 — 10z/ kpc) kms" 1 . The proper motion can be 
written down from the heliocentric velocities from equation 
|T} subtracting the reflex moti ons of the solar Local Stan- 
dard of Rest (LSR) values from lSchonrich et al.l (|2010h and 
assuming a total solar azimuthal velocity according to the 
IAU recommendations at u^.q = 232 kms -1 . Using this the 
simple Galactic streaming terms are: 



Hi = (fixi + aiVg.Q cos I cos b) /(sk) + omasyr 1 
Mi> = (fbXb + CLbV g ,Q sin b sin I) /(sk) + ctmasyr -1 



with 



(3) 



(4) 



Xl = —U sin I + V cos I 
Xt = sin b(U cos I — V sin I) — Wq cos b 

where k ~ 4.74 kms _1 /( pc mas yr _1 ) does the unit conver- 
sion, U and V are the heliocentric velocities from streaming 
minus the assumed motion of the Sun and fi,ai,ci, fb,a,b,Cb 
are the fit parameters. The terms connected to ai give a 
crude estimate for the effects from halo streaming. 

The derived values for DR8 are presented in Table [2] 
The old, erroneous DR8 and DR7 produce similar results. 
Slicing the sample in colour I could not detect any significant 
changes, just the bluemost 10% of the quasars appear to be 
different. Still, cutting for (g — r) > does not evoke notable 
consequences 10%, as well as removing the 10% of objects 
with the largest proper motions does not have any impact. 

The coefficients fi and a; may be understood roughly as 
contamination fractions, but not on accurate terms because 
the distance dependence of proper motions effects a strong 
dependence of the fractions on the assumed distance of stel- 
lar contaminants. The sign in the expected trend s is not 
predictable a priori, as the ISchneider et all l|2010t ) sample 
may harbour itself some residual stellar contamination that 
competes with possible astrometry problems. In both cases 
the main rotation term /, is significantly negative hinting to 
some more misclassified stars affecting astrometry as com- 
pared to the quasar sample. On the other hand the coeffi- 
cients ai that resemble terms expected for the Galactic halo 
are slightly positive, which might be interpreted as a tiny 
surplus of halo stars hiding in the quasar sample. I rather 
understand it as an omitted variable bias: Neither do I know 
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the exact functional shapes and distance distributions of the 
contamination, nor do I have precise terms at hand to cope 
with the possible chromatic aberrations. Exact terms for 
contamination would demand a complete modelling of the 
entire measurement process including a reasonable model to 
the Milky Way, other galaxies and the exact Sloan selection 
function. If I fit alone either di or fi, both variables turn 
out negative and still their counterparts are negative when I 
fit them with this naively derived value. So one could argue 
that fi is getting excessively negative being balanced by a 
slightly positive a» to cope with these uncertainties. 

In the end I cannot achieve more than a rather phe- 
nomenological correction without being able to determine 
the true cause for the observed trends. The naive approxi- 
mation still fulfils its purpose to parametrise the influence on 
Galactic rotation measurements. Throughout the following 
paper I will use the terms for (g — r) > to correct stelar 
proper motions independent of colour, since no significant 
colour dependence was detected. 

For proper motion errors, which influence the distance 
determinations, I use the pipeline values and checked that 
using the values from iDong et al.l l|201lf ) for DR7 instead 
does not significantly affect the presented analysis. 



3.3 Distances 

In this work I make u se of the statistical dist ance correction 
method developed by ISchonrich et all (|201lT ) . The approach 
relies on the velocity correlations over Galactic angles that 
are induced by systematic distance errors and entails a bet- 
ter precision and robustness than classical methods. The 
method can find the average distance in any sample of a cou- 
ple of hundred stars. The only thing required as first input 
is a smooth distance or respectively absolute magnitude cal- 
ibration that has a roughly correct shape in parameters like 
colours and metallicity. The shape has some importance as 
erroneous assumptions would result in an increased distance 
scatter on the sample bins and thus in residual systematic 
trends for different subpopulations. The use of statistical 
distance corrections is mandatory as there are in all cases ex- 
pected offsets by reddening uncertainties, age uncertainties, 
the uncertain binary fraction, mild differences by evolution- 
ary differences on the main sequence, systematic metallicity 
offsets and by the helium enhancement problem, where the 
isochrones appear to be too faint at a fixed colour. 



3.3.1 First guess distances 

Primary stellar distan ces are derived f r om a dense grid 
of BASTI isochrones (|Pietrinferni et all |2004 l2006h that 
was kindly computed by S. Cassisi for t h e pur pose of our 
age determinations in Casagrandc et all (|201 J ). I use the 
11.6 Gyr isochrones which seems an appropriate compro- 
mise bet ween the suspected age of the first disc populations 



(see e.g. lAumer fc Binnevll2009l ; ISchonrich fc Binnevl [20091 ; 
iBensbv et alB oO^ and the age of halo stars. For metal-poor 
stars with [Fe/H] < —0.6 I account for alpha enhancement 
by raising the effec tive [Fe/H] by 0.2 d e x as suggested by 
IChieffi et al.1 l|l99lT) o r IChabover et all (I1992T). In analogy 
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Figure 3. Upper panel: Assumed absolute r-band magnitudes R 
from isochrone interpolation in the metallicity-colour plane after 
statistical distance correction. Lower panel: Distance correction 
factors s/so against the a priori distances sq in the metallicity- 
colour plane with the smoothing from equation |J7}. 



(2008) I let the alpha enhancement go linearly to zero to- 
wards solar metallicity: 



C [Fe/H] + 0.2 
[Me/H] = { [Fe/H] - |[Fe/H] 
I [Fe/H] 



for [Fc/H] < -0.6 

for -0.6 < [Fe/H] < (5) 

for < [Fe/H] 



to 



the observations of lBensbv et all "( 2007 ): Melendez et al.l 



Some of the stars are located bluewards of the turn-off 
point at this age. To cope with these objects, I search the 
point at which the old isochrone is 0.5 mag brighter than 
the corresponding 50Myr isochrone and extrapolate blue- 
wards parallel to the main sequence defined by the 50Myr 
isochrones. I also experimented with using the turnoff point, 
i.e. the blue-most point of the old isochrone, but this in- 
troduces larger shot noise in distances tha t would hamper 
the di stance corrections via the method of ISchonrich et all 
|201ll ). Comparing the expected r-band absolute magnitude 
of each star to the measured reddening corrected ro, I then 
infer its expected distance modulus and hence distance so. 



3.3.2 Distance corrections 

To derive a correction field, I bin the sample in metallicity 
and (g — i)o colour. I also attempted binning the sample 
simultaneously in distance as well, but found no additional 
benefits, probably because the bulk of the sample is further 
away than 0.5 kpc. I also excluded high reddening regions so 
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that variations of reddening with distance are comparatively 
small. 

There is a trade-off between resolution of the smallest 
possible structures in the correction field and the introduc- 
tion of Poisson noise by lower effective numbers of stars con- 
tributing to the estimate at a certain point in the observed 
space. The strategy is to select every 30th star in the sample 
and estimate the distance shift for its 700 closest neighbours 
(I checked that the bin size has no significant effect) using 
an euclidean metric on the parameter space: 



n = 9 A 

u param (g- 



Fe/H]- 



(6) 



For non-halo or respectively metal-rich stars, which offer 
less accurate statistics by their smaller heliocentric veloci- 
ties, this results in some significant scatter, which is reduced 
by smoothing the distance corrections via a weighted loga- 
rithmic mean: 



so exp 



J2i ln i 1 + X est,*) eX P (~ J2 
7 A 5 




(7) 



where so is the a priori distance estimate from isochrone 
interpolation. The first sum runs over all stellar param- 
eter sets i around which the best-fit distance correction 
x = 1 — 1/(1 + /) or respectively the distance mis-estimate 
factor / has been evaluated. The second index j runs over 
used parameters (here metallicity and colour), to evaluate 
the parameter differences Aij between the star in question 
and the mean values of all evaluation subsets and smooth 
via the exponential kernel with the smoothing lengths <jj. 
1 chose f[p e /jjl = 0.1 dex and a^ g _ i - j0 = 0.04 mag. When 
binning the sample a second time for measuring Galactic 
parameters, I fit the distance for a second time to diminish 
the danger of remaining systematic biases. Anyway there is 
no evident trend, e.g. with distance, in the second step and 
the remaining correc t ions a re on the noise level. 

ISchonrich et al.l (|201ll ) showed that in the presence of 
significant intrinsic distance scatter the method a slight 
mean distance underestimate occurs, which is caused by 
stars with distance overestimates populating the edges of 
the fitting baseline and hence acquiring larger weight. This 
bias was found at 5f ~ 0.5cr^, where / denotes the relative 
distance error and cPj its variance. I confirmed the validity of 
the prefactor 0.5 for the sample in use by applying a range of 
Gaussian broadenings to the first guess isochrone distances. 
Though the prefactor is known, it is not straight forward 
to estimate the values of 07 dependent on metallicity and 
colour. On the red side the stars are relatively firmly placed 
on the main sequence, the reddening vector runs quite par- 
allel to the main sequence, and hence the colour and metal- 
licity errors dominate justifying a moderate error estimate 
07 ~ 15%. On the blue end the placement of the stars on 
the turn-off becomes increasingly uncertain and especially 
on the metal-rich side one must expect significant spread in 
ages and hence absolute magnitudes as well. I hence adopt 
a mild increase in distance depending on colour 



— = 0.01 

s 



+0.06 ((g - i)„ - 0.8) 2 (l + i ([Fe/H] + 2.3) ; 



,(8) 



which is larger on the blue side and increases mildly 
with metallicity. I limit the relative distance increase to a 



maximum of 4.5%, which corresponds to a distance scatter 
07 = 30%. 

Fig. shows the finally adopted absolute magnitudes 
(upper panel) and the correction field for the isochrone dis- 
tances (lower panel) in the colour-metallicity plane. Along 
the blue metal-poor stars I require moderate distance 
stretching. This is the turn-off region, where the 11.6 Gyr 
isochrones might be a little bit too young and get shifted 
upwards in addition to the known suspicion of isochrones 
being too faint. On the other hand a bit of contamination 
by subgiants is more likely because of the reduced difference 
in surface gravity and weaker absorption lines. The problems 
with gravities for the most metal-poor stars are indicated by 
the absolute magnitudes bending back, i.e. the most metal- 
poor stars being made brighter again at the same colour. 
Carbon stars constitute another bias having stronger ab- 
sorption in the g-band, making them too "red" in the colour 
assignment, which brings forth a brightness underestimate. 

In principle the method explores the structure around 
the turn-off and the dominant ages among these popula- 
tions, but owing to the scope of this work I defer this to 
a later paper. At higher metallicities the turn-off line ev- 
idently bends upwards, i.e. to redder colours as expected 
from stellar models. On the other hand there is evidence for 
a younger metal-rich population at blue colours, indicated 
by the absence of distance under-estimates contrasting with 
their metal-poor antipodes. Some vertically aligned features 
might point to minor systematics in the pipeline metallicities 
(e.g. under-estimating a star's metallicity makes us underes- 
timate its distance). The little spike near the turn-off around 
[Fe/H] ~ —1.0 could be an unlucky statistical fluctuation. 
If it was real and not caused by a problem in the pipeline 
it would indicate a rather sharp transition in ages with rel- 
atively young stars at slightly higher metallicities opposing 
the very old population in the halo regime. Throughout the 
range of G and K dwarfs there is a quite constant need for 
mildly larger distances than envisaged by naive use of the 
stellar models. 

To test my app r oach for consistency, I experimented 
with the llvezic et al.l (|2008T ) (A7) relation as well and found 
that as expected the distance corrections robustly produce 
an almost identical final outcome despite the highly different 
input. 

3.3.3 Radial velocities and vertical motion of the Sun 

Despite all efforts put into the SEGUE parameter pipeline, 
it cannot be granted that we can trust the radial veloci- 
ties to be free of any systematic errors of order lkms -1 . 
To shed light on this problem, let us examine the verti- 
cal velocity component. A good estimate of vertical mo- 
tion can be derived from the Galactic polar cones, where 
proper motions and hence distance problems have almost 
no impact on the measurement. Luckily SEGUE has both 
the Galactic North and South pole and I use | sin6| > 0.8, 
while checking that a stricter cut at 0.9 does not signifi- 
cantly alter the findings. These polar cones display a sur- 
prisingly low W Q = (4.96 ±0.44)kms _1 . This is 2kms _1 
off the old value from Hipparcos and the GCS. While this 
may of course trace back to some unnoticed stream in the 
solar neighbourhood, it seems appropriate to check for con- 
sistency. 14272 stars towards the North Pole show Wq = 
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Figure 4. Systematic streaming in the SEGUE sample. The left 
hand panel shows the heliocentric U velocity, while the right hand 
panel gives the galactocentric U g , i.e. the real radial motion at 
the position of the star assuming a solar galactocentric radius 
Ro = 8.2 kpc and V 9i q set accordingly to the proper motion of 
Sgr A*. The sample was binned in separate boxes a third of a kpc 
wide, plotted when they contain more than 50 stars. The sample 
thins out towards the remote stars, which causes the increased 
scatter seen particularly in the galactocentric velocities. For both 
panels I restrict the sample to stars with [Fe/H] > — f. 



(4.26±0.52) kms" 1 , while 3482 stars towards the South Pole 
give Wq = (7.84 ±0.85) kms" 1 . Allowance for a variation in 
the line-of-sight velocities: 



W — Svios sin b — Wq 



(9) 



estimates 8vi os = (2.04 ± 0.63) kms" 1 and Wq = (-6.07 ± 
0.56) kms" 1 . The tighter cut in latitude gives the same line- 
of-sight velocity deviation, but Wq = —6.9 kms" . Binning 
in metallicity brings me to the limits of statistical signifi- 
cance, but reveals a rather robust uptrend in the estimated 
vios offset from near to negligible at solar metallicity to of 
order 5 kms" 1 n ear \Fe/H] ~ —2. It i s tempting to argue in 
accordance with IWidrow et ail (|2012T ) that all those trends 
are reaction of the Galactic disc to some satellite interac- 
tion, but such reasoning is contradicted by more of the sig- 
nal coming from the halo stars, which cannot participate in 
such a response. An explanation by streams conflicts with 
the broad metallicity range concerned that would require 
a conspiracy of nature. A simpler explanation might be a 
warp in the disc, but again this conflicts with the metal- 
licity trend. The easiest explanation for the deviation is a 
residual error in the pipeline. To avail this and stay as close 
as possible to the pipeline values, I use henceforth a min- 
imal correction that keeps the discussed trend just below 
significance: 



uios = «io S ,p ip + (-0.2 + 3.3[Fe/H]/2.5) kms" 1 



(10) 



allowing the corrections to vary between —0.2 and 
—3.5 kms" 1 . 



4 DETECTING GLOBAL ROTATION IN 
SEGUE 

Fig.[4]shows the systematic motion for the sample projected 
into and binned in the Galactic plane. At each bin I evalu- 
ate the mean motion in the heliocentric frame (left panel) 
and the mean motion in the Galactocentric frame (right). 
While I would expect some minor systematic motion even 
in the Galactocentric frame by spiral structure and maybe 
influence by the bar, such contributions will be relatively 
small as most of this sample is at high altitudes, where the 



stars with their large random velocities experience less im- 
portant changes of motion by the relatively shallow potential 
troughs of the spiral pattern. Some of the residual structure 
will likely be blurred out by Poisson noise and distance un- 
certainties, so that it is no surprise not to see any apprecia- 
ble signal. As a good sign there are no direction-dependent 
shifts, which would be expected under systematic distance 
errors, which make especially the azimuthal velocity com- 
ponent cross over into the radial velocity determinations. 
In the left panel, however, we see the very prominent rota- 
tional streaming of the Galactic disc, which is central to the 
following investigations. 

4.1 The Standard of Rest revisited 

The strong influence of rotational streaming on the observed 
radial motion makes any determination unreliable that does 
not explicitly take it into account. Before I turn to its pos- 
itive application, I lay out a simple procedure to deal with 
this effect as a contaminant and measure the solar radial 
LSR velocity. In the literature varying offsets of the solar ra- 
dial motion of up to more than 10 kms" 1 from the value de- 
rived from local samples have been claimed. Besides possible 
velocity cross-overs, e.g. the motion in V t ranslating into a 



vertic al motion (for an example see Fig. 7 in Schonrich ct al 



2011 ) invoked by systematic distance errors ( Scho nrich et al 



20111 ). neglected sample rotation can be the reason. In any 
sample that is not demonstrably rotation-free and is not 
symmetric in Galactic longitude the streaming motion off- 
sets the mean heliocentric radial motion and hence the in- 
ferred Solar motion. The error grows with the remoteness of 
the observed stars. 

Fig. [S] shows the heliocentric U velocities for all stars 
(red crosses) from DR8 and stars with [Fe/H] > —0.5 (light 
blue points) against sin a. I set Ro = 8.2 kpc, and the value 
for V g ,Q accordingly to 248.5 kms" 1 . The metal-rich sub- 
sample displays a nearly linear relationship of the mean 
U velocity versus sin a as expected for rotational stream- 
ing from equation fTJ. This streaming motion is far larger 
than the expected offset caused by the s olar LSR motion 
Uq ~ 11.1 kms" 1 (|Schonrich et al.ll2010t ). Even the mild 
asymmetry in a distorts a naive estimate of Uq in disre- 
gard of rotation to Uq — (3.85 ± 0.43) kms" 1 for metal-rich 
objects with [Fe/H] > -0.5 and Uq = (6.77 ± 0.44) kms" 1 
for the entire sample. Consistent with the rotational bias 
the deviation from the standard LSR value is significantly 
larger for the disc sample. To rid my results from this bias, 
I use the estimator 



U(a) 



Uq + 5 l 



(11) 



where the free parameters 9 and —Uq are the effective 
rotation and the reflex of the solar motion, 5u describes 
noise from both velocity dispersion and observational scat- 
ter that should have zero mean independent of a. I obtain 
(9, Uq) = (208.2T3.6, 12.90T0.42) kms" 1 for the metal-rich 
subsample and {9, Uq) = (149.3 ± 3.7, 13.04 ± 0.46) kms" 1 
for all stars. The formal errors are mildly underestimated 
due to the non-Gaussianity of the U velo city distribution 
and a non- weighted fit. This differs from the lSchonrich et al.l 
(|2010T ) estimate at marginal significance. 

Yet there are larger uncontrolled systematic biases: This 
fit implies that the effective rotation term carries away all 
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Figure 5. U velocities versus sin a for all stars in the SEGUE 
sample (red crosses) and objects at disc mctallicities ([Fe/H] > 
—0.5, blue points). The line represents the linear fit for the metal- 
rich subsample, the dark green error bars indicate the mean ve- 
locities in bins of sin a. The lower panel is a close-up of the central 
region. 



the observed bias on U®. While the detailed rotation pat- 
tern is not of interest here, the measurement will still be 
biased if one encounters a larger fraction of slow rotators on 
one side than on the other. E.g. a different number of halo 
stars or unmatched differences of rotation with galactocen- 
tric radius can distort 8, which correlates with U®. In this 
sample U® shrinks by 0.5 km s for every 10 km s -1 decrease 
in 0. In the present case such effects are below detectability, 
as can be seen from the lower panel in Fig. O where the 
binned means (green error bars) show no significant devi- 
ation from the linear relationship. Further, any remaining 
systematic distance errors can lead to systematic velocity 
shifts. The first error can be controlled by selecting only 
stars with small a, the second by selecting stars near the 
Galactic centre and anticentre directions. In principle this 
could be tested by geometric cuts on the sample, however, 
there are no data near the Galactic plane. Exploring the 
full parameter space of possible cuts reveals some expected 
fluctuations of order 1.5 kms -1 , but no striking systematics. 
For example |sinaj < 0.1 gives U® = (13.29 ± 0.54) km s -1 
in concordance with the full sample. 

More importantly, cuts for a smaller distance range, 



which is more reliable concerning the proper motion un- 
certainties favour U® ~ 14 kms -1 . The difference is caused 
by the ~ 15% most remote stars displaying a rather ex- 
treme U® ~ 9 kms -1 . E.g. cutting for s < 2.5 kpc I get 
U® — (14.26 ± 0.50) kms -1 , again consistent with U® = 
(14.08 ± 0.57) kms -1 for |sina| < 0.1, demonstrating that 
this drift in distance is not an effect of some rotation error. 
A similar rise in U® to up to 15 km s -1 can be achieved by in- 
creasingly tight cuts on the signal to noise ratio, which indi- 
rectly selects a shorter mean distance. A global line-of-sight 
velocity increment 5v\ OB — 2 kms -1 increases the estimate 
for U® by about 0.5 kms - . 

It is interesting to investigate into the question of a 
trend in ra dial velocities with galactocentric radius R as 
suggested in ISiebert et ai]||201lf ). To this purpose I estimate 
the equation 

U = 6 sin a + r)(R - Rq) - U® + e v (12) 

where 9 is the effective rotation, R the galactocentric radius, 
77 = d{U)(R)/dR is the effective (heliocentric) radial veloc- 
ity gradient on the sample, and ev is the noise term from 
velocity dispersion. For convenience I limit sin a < 0.1 
as before and get on the full distance range rj = (1.53 ± 
0.55) kms -1 kpc -1 at U® = (14.22 ± 0.64) kms -1 . On the 
shorter distance range of s < 2.5 I get U® = (14.22 ± 
0.67)kms -1 and 77 = (0.29 ± 0.7) kms -1 kpc -1 , deliver- 
ing no sensible gradient at all. Restricting [Fe/H] > —0.5 
I get r) = (1.05 ± 0.7) kms -1 kpc -1 . This residual trend 
is caused by the 1000 stars with [Fe/H] > measuring 
rj = (3.7 ± 1.9) kms - kpc -1 . While this may finally be the 
desired signal, it is puzzling to have no more metal-poor 
counterpart. Hence this is more likely either be a statisti- 
cal fluctuation or some flaw in the sample. Parameters and 
distances for stars with super-solar metallicity in SEGUE 
are likely problematic, particularly as I do not have suffi- 
cient numbers to constrain their mean distance very well. 
It has to be noted that this sample is at rather high alti- 
tudes and is dominantly in the outer disc, so that velocity 
trends in objects of the inner thin disc might not show up. 
Hence the missing detection of an appreciable trend (and ex- 
clusion of a significant devia tion from circularity) does not 
disprove lSiebert et alj (|201ll ). More of a problem for their re- 
sult should be the larger value of U® , which likely reduces or 
diminishes their estimated gradient, as their analysis mainly 
relies on assumption of the old value and does not have suf- 
ficient data points in the outer disc to provide a decent fix. 

After all, when accounting for the systematic uncertain- 
ties in the data, I cannot detect any reliable hint for a mo- 
tion of the solar neighbourhood against some more global 
Galactic Standard of Rest. Combining these estimates with 
the result from the mean direction of motion (see Section 
I4,5[l of (13.9 ± 0.3) kms -1 with a systematic uncertainty of 
~ 1.5 kms -1 and the reader may either use this value or 
a lower U® = (13.0 ± 1.5) kms -1 , taking into account the 
value to the older GCS. 

4.2 Divide and conquer: The rotation of 
components 

So far I exclusively looked at the global rotation pattern 
and its influence on U® , but the new rotation estimator can 
do more: Any measurement relying on observed azimuthal 
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Figure 6. Measuring the intrinsic rotation of the populations. 
I assume Rq = 8.2kpc and V 9 ,0 = 248.5kms _1 . After sorting 
the sample in [Fe/H] I move a ~ 2500 stars wide mask in steps of 
1250 so that every second data point is independent. The red data 
points show the measured rotation speed in each subsample, while 
the blue data points with small errors give the mean azimuthal ve- 
locity (V g ) in each bin. The green line shows (V g ) w = . wiV g ^ 
from equation H13I I. which is weighted with (sina — sinct) 2 . 



velocities does not provide the real rotation of the compo- 
nent in question, but only a value relative to the Sun, from 
which we conclude its azimuthal speed by discounting for 
the Galactic circular velocity plus the solar azimuthal LSR 
motion (plus a small part of the solar radial motion when 
\a\ > 0). On the contrary the radial velocity part of (QJ uses 
the velocity differences at a range of Galactic positions and 
gives the "real" rotation of the component without signifi- 
cant influence by assumptions on the motion of the Sun. The 
measurement is in most cases less precise than the azimuthal 
measurement, but is nevertheless interesting by its different 
set of assumptions. Consider a possible presence of stellar 
streams: They can distort azimuthal velocity measurements 
soliciting an independent source of information: In general 
a stream will then show up in highly different values for the 
two rotational velocity indicators. The price one has to pay 
for the radial velocity estimator is a larger formal error, the 
need for extended samples, a different dependence on the 
assumed Ro (which is usually a bonus) and some stronger 
vulnerability to distortions from axisymmetry, e.g. influence 
by the Galactic bar. 

As discussed by ISchonrich et al.l (|201l|), dissecting a 
sample by kinematics is detrimental for this kind of kine- 
matic analysis, but one can divide the Galaxy into its com- 
ponents by a selection in metallicity. 

Fig. [6] shows the rotation measurement from the ra- 
dial velocity estimator when I dissect the sample into bins 
of ~ 2500 objects in metallicity and move half a bin in 
each step (large red error bars). The two plateaus of disc 
and halo rotation with the transition region in between 
show up nicely. The plateau of disc rotation stays about 
15 kms -1 below the circular speed of the disc due to asym- 
metric drift. Clearly there is no significant average net rota- 
tion in the metal-poor Galactic halo. Despite the far smaller 
formal error, this would be more difficult to state on firm 
grounds by using the average azimuthal velocity alone (blue 
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Figure 7. Demonstrating the detectable effect of errors in Galac- 
tic parameters. Settings are identical to Fig. HI] with one difference 
each. In the top panel I set V 9} q = 228.5 km s -1 close to the IAU 
standard value. While the measured rotation speed from U stays 
identical, the determination by V g is shifted down into stark dis- 
agreement. In the bottom panel I set Ro = 7.0 kpc and lower the 
Solar motion to F 9 ,0 = 225.5kms -1 . The compression in the 
measurement from U velocities compared to the V g measurement 
is clearly visible. 



error bars), which does depend strongly on the solar Galac- 
tocentric motion. In fact the azimuthal velocity estimator 
has been shifted onto the U estimator's curve by matching 
V 9 ,0 = 248.5 kms -1 , which turns us to the topic of the next 
subsection. 



4.3 Deriving Galactic parameters 

More important than mutual confirmation of results one can 
play off the two rotation measurements (from U and V g ) 
against each other. As stated before, only V g depends on 
the solar azimuthal velocity, so that forcing both rotation 
estimators into agreement we obtain an estimate for the to- 
tal solar azimuthal velocity in the Galactic rest frame V^q. 
By subtraction of the solar LSR motion this directly esti- 
mates the Galactic circular speed V c near the solar radius 
without any need to invoke complicated models. The appar- 
ent difference caused by a wrong Vg,e is demonstrated in 
the top panel of Fig. [7] For this plot I lowered V 9t @ from 
the optimal value by about 20 kms -1 . Comparison with the 
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previous figure reveals that the red bars are unchanged as 
expected, while the azimuthal velocity based measurements 
are pushed upwards into stark disagreement. 

To bring the two rotation measurements on common 
ground one has to apply weights to the estimate of the mean 
azimuthal velocity (V g ) for each subpopulation: In a least 
squares estimate on the rotation speed 9 each star appropri- 
ates a weight proportional to its distance from the popula- 
tion mean on the baseline d a — | sin a — (sin a) j . As the value 
Uh invoked by rotation varies linearly on this baseline, each 
star attains the weight w — (sin a — (sin a)) 2 when measur- 
ing rotation. Introducing a simultaneously weighted mean 
for the mean azimuthal velocity 



(v g ) w =J2 wiV »> 



(13) 



where the sum runs over all stars i, any structural differ- 
ences between the two estimators are eliminated, since every 
population enters both estimators with equal weights. This 
measurement of the Solar velocity does thus not require any 
modelling of Galactic populations or any other deeper un- 
derstanding at all. The procedure is then very simple: 



• Assume -Ro and a first guess solar velocity vector 



U, 



• Estimate 9 and (V g ) 

• R( 

9, (V g ) w until convergence 



Reduce Vg s by the difference (V g ) w — 9 and recalculate 



The only remaining uncertainty (apart from systematic 
measurement errors) can arise if the geometry of the sam- 
ple allows for a distortion in the radial velocity estimator. 
I hence recommend controlling the value of Uq, which is 
well-behaved for this sample, and in other case fixing it to a 
predetermined value. In the latter case the weights of (V g ) 
must be changed to w — (sin a) 2 . 

The weighted averages on the azimuthal velocities are 
drawn with the green line in Fig. [6] Generally the weighted 
values in the disc regime are a bit lower, because the ex- 
treme ends on sin a are populated by remote stars at higher 
altitudes and at larger asymmetric drift outweighing other 
effects like the polewards lines of sight reaching higher alti- 
tudes. 

Via the Galactic angle a the assumed solar galactocen- 
tric radius enters the derivation. As depicted in Fig. [2] the 
larger Ro, the smaller the derived |q| will be for stars on our 
side of the Galaxy. This results in a larger rotation speed 
9 required for the same observed pattern and so finally in- 
creases the estimate of V g ,@. Fortunately the rotation speed 
of Galactic components is lower than the solar value and 
further - in contrast to the linear Oort constants - the de- 
pendence on a is not linear in distance, which helps towards 
a relationship natter than a mere proportionality between 
V c and Vq. This is shown with the red error bars in Fig. [8] 
For the practical calculation I bin the sample in distance and 
move a mask of ~ 2500 stars width in steps of 625 stars. This 
stabilises the results against scatter from the coarse binning, 
the error is corrected for the induced autocorrelation. The 
result is then drawn from a weighted average of the distance 
bins. The binning in distance is mandatory in this sample 
not mainly because of reddening uncertainties, but because 
of the previously described anomaly in the mean radial ve- 
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Figure 8. Measuring the total azimuthal velocity V 9i q =V c + Vq 
of the Sun by forcing an agreement between the radial velocity 
based rotation estimator and the azimuthal velocities. For each 
assumed Ro the sample was binned in distance at bins of ~ 2500 
stars and then a weighted average taken, which is shown with red 
error bars. To guide the eye, I add the lines of a linear fit. For 
comparison I plot the constraint from the proper motion of Sgr 
A* with blue lines. 
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Figure 9. Using the trend in determined rotation speed versus 
mean azimuthal velocity to get another constraint to the galacto- 
centric radius. For this purpose I sorted the sample by metallicity 
into bins of 2500 stars and estimated dVg t Q/dV g (vertical axis in 
the plot, a dimensionless quantity) while varying Ro- 



locities with distance. This way contaminations are kept at 
a minimum level. 

By using the proper motion of Sgr A* (deter mined as 
6.379±0.026masyr _1 bv lReid fc Brunthalerll20oi ) as a sec- 
ond constraint the Galactocentric radius can be fixed to 
Ro = 8.26^;^ kpc and V 9i0 = 250±iJ kms~\ For the re- 
sulting circular speed one has to subtract about 12kms _1 . 

One can even gain an estimate for the Galactocentric 
radius that is fully independent of any other results. For 
this purpose we have to recapitulate again the impact of 
Ro on the derived Vg,© for different populations, shown in 
Fig. [2] The measurement of V S: q compares the weighted az- 
imuthal velocity in any population to the derived absolute 
rotation 9. Changing the Galactocentric radius Ro makes 
the Galactic coordinate system bend differently and hence 
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gives only minor changes to V g , but strongly affects sin a. 
This way 9 is nearly proportional to Ro. A weakly rotat- 
ing population has a far smaller value of 9 and hence ex- 
periences a far smaller change than fast rotating disc stars. 
If we now underestimate -Ro, the estimates for 9 are com- 
pressed compared to the V g based estimator as demonstrated 
in the bottom panel of Fig. [7] In other words the estimates 
of V 9t Q for disc stars get sheared against the estimates for 
halo stars. However, for all populations one measures the 
same quantity, i.e. the azimuthal speed of the Sun, which 
cannot change between the samples. A deviation in Ro is 
hence easily detected by dissecting the sample in metallic- 
ity, which is (as can be seen from Figure [6} a very good 
proxy for rotation. I cut the sample again into bins of ~ 2500 
stars - now in metallicity instead of distance - and evaluate 
the trend V 9 ,q dV gt Q/d(V g ) w of V 9 ,q against the mean az- 
imuthal speed between the bins by fitting the linear equation 
Vg,© = (dV g ,q / d{V g ) w ) (V g ) + Vo with Vo as additional free 
parameter. From the red error bars in Fig.|9]one can see that 
already on the quite modest sample size and extension of the 
~ 50000 usable stars in DR8 I get a useful, independent esti- 
mate of the solar Galactocentric radius Ro = 8.29^'^ kpc. 

As in the first measurement I used the absolute value 
of rotation, while here the unit-less slope dV g$ @/d{V g ) w is 
employed, the two results are formally independent and 
can be combined to Rp = (8.27 ± 0.29) kpc in excel- 
lent agreement with iMcMillanl (|201lh . This translates to 
y q Q = (250 ± 9) km s -1 . Working in the LSR value of 
ISchonrich et al.1 (|201Ch I hence obtain V c = (238±9) kms" 1 . 



4.4 Assessment of systematic errors 

To assess the systematic uncertainties I vary some of the 
assumptions. First I checked that the choice of bin size does 
not change the results more than the general noise induced 
by bin changes as long as the bins are not unreasonably 
small. Galactic rotation at fixed Ro is increased by about 
2 kms -1 , when I globally augment distan ces by 10% and 
would be about 4 kms -1 lower without the lSchonrich et al.l 
(|201lf ) distance correction. Increasing the distance scatter by 
applying an additional 10% raises V g ,@ by 0.6 kms -1 at fixed 
radius, application of a 30% error lowers it by 1.4 kms -1 . 
The mild lowering is expected, since the strong distance 
scatter should give a global distance underestimate of 4.5% 
(consistently the rotation measurement from U velocities 
drops by ~ 2 kms -1 ). It is obvious that such a scatter has a 
minor effect, as the relation between Uh and sin a is linear 
and hence at small changes the mean of Uh is only mildly 
affected. This is a benefit from operating remote from the 
Galactic centre where the changes in Galactic angle become 
more rapid. The mean azimuthal velocity is mostly affected 
by distance scatter via the bias on average distance derived 
from our statistical method. Tests reveal an expected sys- 
tematic bias of order lkms -1 . 

Similarly the slope of the derived rotation speed 
dV 9t Q /dVg depends mostly on distance. Here the dependence 
is larger: It reacts by 0.08, i.e. about 1.5 standard deviations 
to a forced increase of 10% in distance and by —0.08 to a 
forced decrease of 10%. Removing the statistical distance 
correction would lift values by about 0.2. Increased distance 
scatter without distance correction increases the value by 



0.011 for an additional scatter of 10% and lowers it by 0.023 
for a scatter of 30%. 

Gravity cuts have minor influence on my results: The 
mean rotation rate is lowered by 0.9 kms -1 when I lift the 
minimum gravity to log(g) > 4.1 and at the same time the 
trend in azimuthal velocity rises by 0.02. Tighter cuts are 
not feasible for the trend, as I lose all metal poor stars. Re- 
placing the sloping cut by a fixed log(g) > 4.0, did not signif- 
icantly alter the results ether. While the distance estimator 
strongly hints to some contamination of the dwarf sample, 
the observed deviation in Ro is still not significant and likely 
results partly from random scatter, as many metal-poor ob- 
jects are lost by this tighter cut. The robustness of the result 
against changes in the gravity selection points to the statis- 
tical distance estimator being able to cope with the different 
contamination by adjusting the mean distances. 

The largest uncertainty appears to arise from the proper 
motion determinations. The correction I apply raises rota- 
tional velocities by about 6 kms -1 and in light of the unsat- 
isfactory physical reasoning as discussed in Section f3. 21 it is 
uncertain how well the correction works. A little reassurance 
is given by the second estimator, its unitless slope rises by 
only 0.01 upon removal of the proper motion correction. The 
applied line-of-sight velocity correction lowers V 3i q by about 
1.2 kms -1 and lowers the unitless slope by just —0.008. The 
better robustness of the second estimator for Ro reflects the 
fact that any spurious rotation affects both halo and disc 
rotation in a similar way and hence their difference is more 
robust than the mean. 

Further I find that a change of Uq in the calculation 
of galactocentric azimuthal velocities (a minor bias is ex- 
pected in the lopsided sample) is irrelevant to the estimated 
quantities. However, if the rotation estimates themselves are 
done with fixed Uq instead of using it as a free parameter, 
the estimates agree fully when setting Uq — 13.9 kms -1 , 
and V 3> q falls by about 4 kms -1 for every lkms -1 reduc- 
tion in Uq. Depending on how much integrity we ascribe to 
the astrometry in the sample this sets an interesting rela- 
tion between Uq and Ro. Similarly Uq can be determined by 
demanding the estimate for V g> Q to be consistent in Galac- 
tic longitude I. Not surprisingly the favoured value is again 
Uq ~ 14 kms -1 . I further tested that there are no trends of 
the rotation estimate (with fixed Uq) with colour, distance 
or latitude. 

The estimator of Uq from Section 14.11 depends weakly 
on distance. Uq increases by about 0.8 kms -1 for a global 
10% increase in distances and decreases by about the same 
amount for a distance increase. Its counterpart 9 changes by 
around 3 kms -1 with some fluctuations as different distance 
cuts can affect the sample composition. The proper motion 
correction is responsible for about 6 kms -1 of the measured 
9 and 0.4 kms -1 of Uq, while the line-of-sight velocity cor- 
rection acts in the opposite direction, lowering 9 by 2kms - 
and Uq by ~ 0.2 kms -1 . 



4.5 Using the direction of motion 

As a third strategy for fixing Galactic parameters I suggest 
the direction of motion. All that needs to be done is to com- 
pare the expected Galactic angle a to the angle that the 
Galactocentric mean motion in the subsamples has against 
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Figure 10. The upper panel shows the angle a v from equation 
IU4H against the position angle a (in rad) of the mean motion 
vector in the plane for bins with more than 10 stars. The distri- 
bution of bins is irregular because the mean values in x and y 
were used in contrast to Fig. [4] where bin boundaries were used. 
The lower panel shows the resulting estimates for Rq varying V 9j 
and allowing for Uq as free parameter. 

the local azimuth. The upper panel in Fig. \W\ shows the 
angle 

f ( (U) + Ue,o \ m x 
0h ,=arctan ^ v)+Vg[a j (14) 

of the mean velocity for a binned sample of all stars with 
[Fe/H] > —1.0 (above that value I expect sufficient rotation) 
assuming a solar azimuthal velocity V 9 ,q = 248.5 kms -1 , 
Uq,o = 13.0 kms -1 and plotting only bins with more than 10 
stars. The lines show positions of constant Galactic angle a. 
For easier fits I note that the ratio {(U) + Uq,q) /((V) + V 3> q) 
behaves as y/(x + Ro) as long as we stay away from \a\ ~ 
(avoiding the arctangent gives better convergence when 
starting far from the optimal parameters). It is obvious that 
this kind of statistic critically depends on the choice of Uq 
and the total azimuthal velocity of the Sun. Formally the 
agreement between the position angle and the velocity angle 
is optimised by fitting 

a'(Ro,U e ) = arctan (-^-) +j(U Q ) (15) 



7 (tfe) = arctan [^y^ ) ~ arctan [-^y^ ) (16) 

to a v , where ^/(Uq) does the correction from the first guess 
Uofi to the parametrised solar radial velocity Uq and x and 
y are the coordinates of the local Cartesian frame in the 
radial (outwards) and azimuthal direction. 

While I can fit Uq directly to these data, this is not 
possible for Vtot as the fit would then converge to the wrong 
global minimum: Let Vtot and Ro go to infinity and the fit 
becomes perfect with all ratios/angles zeroed. In this light 
I vary Kot externally and allow for Uq and Ro as free fit 
parameters. All bins were weighted by the inverse number 
of stars they contain and I only accepted bins with more 
than 10 objects. I note that velocity dispersions potentially 
exert a (predictable) distortion on a v by its non-linearity, 
but found this effect to be very small on the current sample. 

It is not surprising that the result for the solar radial 
velocity is Uq = (13.84T0.27) kms -1 in excellent agreement 
with the direct measurement from stellar U velocities and 
with Section l4~4l 

The result for Ro is shown in the lower panel of Fig. 1101 
where I plot the new datapoints from angle fitting in addi- 
tion to the previously discussed data. The formal error on 
Ro at fixed V 3i q is a small 0.12 kpc, but their relation largely 
resembles that of Sgr A*. I would need a larger sample ex- 
tension and/or significantly more stars to compete with Sgr 
A*, which will be achieved with Gaia. For now, formal in- 
clusion of the third constraint hence does not change the 
final results and just reduces the formal errors a bit. Never- 
theless the potential of this estimator for larger samples is 
evident. The residuals of the fit will likely outline perturba- 
tions of the Galactic potential and substructure with data 
from Gaia, but in the present sample I could not detect any 
structure. 

Adding stars with halo metallicities to the test just 
drives up the error bars by a factor of 2 because of the higher 
dispersion, while no rotating stars are added. Ro stays the 
same while Uq shifts to (14.10 ± 0.41) kms -1 . Restricting 
the disc sample to a distance closer than 2 kpc increases Ro 
and Uq by insignificant 0.06 kpc and 0.06 kms -1 . The value 
for Uq depends on the position of the Galactic centre. If our 
coordinate system was off-centre by 1 degree in Galactic 
longitude, Uq could shift down by about 3 kms -1 . 



5 CONCLUSIONS 

The most important outcomes of this work are three 
modelling-free and simple estimators for the Solar azimuthal 
velocity and hence the local circular speed V c and for the So- 
lar Galactocentric radius Ro. 

On this course I developed the idea that in a spatially 
extended sample the absolute rotation of stellar components 
can be measured from systematic streaming in the heliocen- 
tric radial direction. The stars on one side of the Galactic 
centre show an opposite heliocentric U velocity to those on 
the other side. This value has a lower formal precision than 
the classically used azimuthal velocity, but can boost accu- 
racy by its relative independence from assumptions about 
the velocity of the Sun. 

The rotation in any extended sample severely affects 
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determinations of the Local Standard of Rest: Uq and Wq 
are frequently determined via simple sample averages in 
each component. The rotation bias in Uq gets more im- 
portant with increasing distance and impacts all presently 
available big surveys, since they are lopsided, i.e. asymmet- 
ric in Galactic longitude. For SEGUE dwarfs it amounts 
to ~ 10 kms -1 . Accounting for rotation the otherwise ob- 
served difference between disc and halo stars disappears 
and combining this classic method with the estimate from 
the mean direction of motion described below, I find Uq = 
(14.0±0.3) kms -1 with an additional systematic uncertainty 
of about 1.5 kms - . The value is ~ 3 kms -1 larger than 
from the Geneva-Copenhagen Survey, but still within the 
error margin. While the GCS is clearly affected by stellar 
streams, the presented value may be distorted by the prob- 
lematic Sloan proper motions and possibly residual distance 
errors. 

While there could of course be some interesting physics 
involved, a systematic difference of ~ 4 kms -1 in the 
average W motion between cones towards the Galac- 
tic North and South Poles points to a systematic er- 
ror in the line-of-sight velocities by ~ 2kms . This is 
not implaus ible in light of the adhoc shift of 7.3 km s -1 
appli ed in |Adelman-McCarthv et ah! l200Sl : lAihara et al.l 
l201ll ). The correction reconciles Wq to reasonable agree- 
ment with Hipparcos and the Geneva-Copenh agen Survey 
|Holmberg et al.1 l2009t lAumer fc Binnevl l2009h albeit at a 
lower value of about 6 kms -1 . 

Comparison of the absolute rotation measure 9 based 
on heliocentric U velocities to the mean azimuthal velocities 
V g in a sample delivers the solar azimuthal velocity V 3i q. 
This measurement is correlated with the assumed Galac- 
tocentric radius Ro. Combining this relation with another 
datum like the proper motion of Sgr A* one can determine 
both Ro and V g ,@. For DR8 I obtain Ro = 8.26±oH kpc 
and V g ,Q = 250^5 kms -1 . By dissecting the sample via 
metallicity into slow and fast rotating subgroups I can in- 
dependently infer the Galactocentric radius from their com- 
parison: A larger Ro reduces a and hence nearly propor- 
tionally increases 9. Thus fast rotators experience a larger 
absolute change in the rotation speed 9, putting their value 
of V g ,Q at odds with that from the slow rotators when Ro is 
wrong. Enforcing consistency provides Ro = 8.29io'«;4 kpc, 
and in combination with the simple rotation measure and 
the proper motion of Sgr A* I get R = (8.27 ± 0.29) kpc 
and V g ,Q = (2 50 ±9) kms -1 in ex cellent agreement wit h 
the values from iMcMillanl (|201lh or lGillessen et al l <|2009h . 
The circular speed V c = V s ,q — Vq using the LSR value o f 
Vq = (12.24 ±0.47 ±2) kms" 1 from ISchonrich et all (|2010T ) 
is then V c = (238 ± 9) kms" 1 . 

The third approach uses just the direction of motion 
in the Galaxy to estimate Ro- Assuming near axisymme- 
try again the direction of motion in a sample should always 
point in the direction of the azimuth. Fitting the Galactic 
angle throughout the plane to the angle of motion provides 
Ro. This measurement displays a strong dependence on the 
solar radial motion apart from the dependence on V 3i q. For- 
tunately Uq and Ro are not degenerate and I obtain again a 
relatively large Uq = (13.84 ± 0.27) kms -1 , confirming the 
previous conventional analysis, but with an even smaller for- 
mal error. As the resulting relationship between V g ,Q and Ro 
is only weakly inclined against the result from Sgr A* , this 



approach will be of relevance for larger samples and here 
just provides some reassurance. 

The laid out strategies are extremely simple. They do 
not rely on any modelling with possible hidden assumptions, 
so that possible biases are readily understood. They rely 
on the sole assumption of near axisymmetry of the Galac- 
tic disc. They are designed for thick disc stars that reside 
high above the Galactic plane and are by their position and 
high random energy far less affected by perturbations of the 
Galactic potential, like the spiral pattern. Their modelling 
independence is supported by the fact that by construction 
none of the estimators depends on the detailed changes of 
sample composition with position in the Galaxy or on the 
radial dependence of V c . Still some uncertainties must be 
pointed out: The methods are vulnerable to large-scale dis- 
tortions from axisymmetry in the Galactic disc. The bar 
region should be avoided, and there may be a residual sig- 
nal from spiral arms. I could not detect significant struc- 
tures connected to this, so the consequences should be rather 
small, particularly as I stay out of the region dominated by 
the bar and as the sample is sufficiently extended not to 
cover just one side of a spiral arm. 

The currently available data add to the possible bi- 
ases by their uncertain distances, radial velocities and espe- 
cially astrometry. The Segue proper motions display a dis- 
tinctive systematic pattern on quasar samples. While I use 
the quasar catalogues for corrections, the solution remains 
unsatisfactory: I could not quantify separately the physi- 
cal reasons, i.e. chromatic aberrations, astrometric "frame- 
dragging" by stars mistaken for Galaxies or even focal plane 
distortions of the telescopes, and also there might be some 
undetected dependence on the stellar colours. Further some 
minor uncertainty in parameters is caused by possible biases 
in line-of-sight velocities. 

I m ade extensive use of the distance corrections devel- 
oped by ISchonrich et "all (|201ll V Besides that this project 
would have been futile without the accuracy achieved by 
the statistical corrections, the outcome is prone to all the 
weaknesses of that method. Especially streams and wrong 
assumptions about the velocity ellipsoid can induce system- 
atic distance errors of order 5%. The contamination in the 
sample will vary to some extent with Galactic coordinates, 
an effect that I did not model. This can lead to a bias, be- 
cause the distribution of weights in the rotation estimators 
and the distance estimator is different. The accuracy could 
be far better could I use more stars and had I better control 
of the systematic errors by another dataset. Part of this will 
be resolved by use of additional samples like RAVE. Fortu- 
nately the estimators show very different behaviour on the 
different biases: The proper motion problem affects almost 
exclusively the first estimator for rotation with a correction 
to V g ,Q of ~ 6 kms - , while the difference between fast and 
slow rotators is nearly untouched. Vice versa, the difference 
estimator reacts strongly to distance changes, while the first 
estimator does not - a 10% change shifts it by ~ 2 kms -1 . 
Despite their benign distribution the systematic uncertain- 
ties on Vg t Q and Ro may come close to the formal errors. 

If the reader should take only one point from this note 
then let it be this: With the advent of the large surveys 
stellar samples are regaining their place as a primary source 
to obtain not only the Local Standard of Rest, but global 
Galactic parameters. Already a sample of ~ 50000 stars from 
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SEGUE provides formal accuracies for the galactocentric ra- 
dius and the solar azimuthal velocity that are competitive 
with any other known approach and without any need for 
modelling. 
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